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Abstract 

Quasi-elastic helium atom scattering measurements have provided 
clear evidence for a two-dimensional free gas of Xe atoms on Pt(lll) 
at low coverages. Increasing the friction due to the surface, a gradual 
change of the shape of the quasi-elastic peak is predicted and analyzed 
for this system in terms of the so-called motional narrowing effect. 
The type of analysis presented here for the quasi-elastic peak should 
be prior to any deconvolution procedure carried out in order to better 
extract information from the process, e.g. diffusion coefficients and 
jump distributions. Moreover, this analysis also provides conditions 
for the free gas regime different than those reported earlier. 



1 Introduction 

One of the main observables in adsorbate diffusion on metal surfaces is 
the full width at half maximum (FWHM) of the quasi-elastic (Q-)peak as a 
function of the parallel momentum transfer, AK = Kj — Kj. Different mag- 
nitudes of interest can be easily extracted from it, e.g. diffusion coefficients 
and jump distributions. In the 1950s van Hove [1] and Vineyard [2] studied 
the broadening/narrowing of the Q-peak for several simple models (ranging 
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from the free particle to a lattice atom and the liquid phase) in order to 
characterize the nature of the phase under study in terms of this property 
(if possible). However, no attention has been usually paid to the actual 
lineshape of the Q-peak. Is the shape of the Q-peak preserved when the 
surface friction and the parallel momentum transfer of the particles probing 
the adsorbates change? This is a critical question; a precise knowledge of the 
natural lineshape is crucial in the experimental deconvolution procedure for 
a better description of the diffusion mechanisms underlying the adsorbate 
dynamics. Nowadays, the wealth of experimental data available in surface 
science (obtained either by standard time-of-flight techniques O d] or by 
spin-echo measurements [5]) allows us to tackle this question successfully, 
as we show in this work. In particular, quasi-elastic helium atom scattering 
(QHAS) has provided [4] the first experimental evidence for a fully mobile 
two-dimensional gas of Xe atoms on Pt(lll) at low coverage {0 = 0.017), low 
incident helium atom energy (Ei = 10.15 meV) and a surface temperature 



One of the theoretical approaches used to interpret the experimental 
results from QHAS measurements is the classical Langevin formalism [6l 
[7], sometimes combined with standard molecular dynamics techniques [1]. 
Here we are going to use such an approach in order to analyze the Xe- 
Pt(lll) system at the experimental conditions considered by Ellis et al. 
[4], analyzing the Q-peak lineshape in the light of the so-called motional 
narrowing effect [6l [H |9]. As will be seen, this type of analysis should be 
prior to any deconvolution procedure aimed to better extract information 
from the process, such as diffusion coefficients or jump distributions. 

The organization of this work is as follows. In order to be self-contained, 
in section [2] we present the theoretical approach that we have followed to 
obtain the analytical formulae that will allow us to interpret our numerical 
results; the numerical methodology is also presented in this section. The 
results derived from our approach are shown and discussed in section [3l 
Finally, the main conclusions are summarized in section HI 

2 Theoretical approach 

From QHAS experiments one obtains the differential reflection coefficient 
which, in analogy to liquids [H [TO], can be expressed as 
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and gives the probability that the He atoms scattered from the diffusing 
collective reach a certain solid angle with an energy exchange fkv = Ef — Ei 
and parallel (to the surface) momentum transfer AK = Kj— Kj. In the right 
hand side (r.h.s.) of equation ([1]), is the adparticle concentration; is 
the atomic form factor, which depends on the interaction potential between 
the probe atoms and the adparticles; and S{A'K,uj) is the so-called dynamic 
structure factor or scattering law, which is the observable magnitude in this 
type of experiments and provides complete information about the dynamics 
and structure of the ensemble of adparticles, and therefore also about the 
surface diffusion process. For example, information about long distance 
correlations is obtained from ^(AK,^) when considering small values of 
AK, while long timescale correlation information is available at small values 
of Huj. 

The starting point of our approach consists of expressing the dynamic 
structure factor as [H El [TO] 

S{AK,uj) = j e-''^^ I{AK,t)dt, (2) 

where 

is the intermediate scattering function, where the brackets denote (ensemble) 
averaging over the adsorbates' trajectories, R(t). In vak is the adpar- 
ticle velocity projected onto the AK-direction (note that AK = ||AK|| is 
the length of AK). 

Both the dynamic structure factor and the intermediate scattering func- 
tion can be readily obtained from Langevin numerical simulations as fol- 
lows. Within a low coverage regime (as happens in Ellis et al. experiment 
[1], where 9 = 0.017), adsorbate-adsorbate interactions can be neglected, 
and diffusion is fully characterized by only studying the dynamics of an iso- 
lated adsorbate on a metal surface. This is the so-called single adsorbate 
approximation. The adsorbate-substrate interaction includes the effects due 
to the surface corrugation as well as those arising from the surface thermal 
vibrational modes. Considering a temperature dependent expansion, both 
effects can be separated. Thus, one the one hand, there is an adiabatic, 
periodic adsorbate-substrate interaction potential V which is temperature 
independent, namely the zero-temperature potential (i.e. Tg = K). On the 
other hand, there is a coupling term accounting for the vibrational effects 
induced by the temperature on the (surface) lattice atoms that act on the 
adsorbate, which can be substituted by a stochastic noise source [IT]. This 
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allows to use a Langevin approach [3] , where the force acting on the adsor- 
bates is given by two contributions: (1) a deterministic force F = — W, and 
(2) a stochastic force G{t). As in standard Langevin molecular dynamics 
simulations OH], the stochastic force has features of a Gaussian white noise 
(the diffusion process is considered as a Brownian-like motion [H]), i.e. 

{G{t)) = 0, (4) 
(G(ti)G(t2)) = 2mjkBTs5{t2-h), (5) 

where m is the adsorbate mass and 7 is the adsorbate-substrate coupling 
strength or friction coefficient. 

Taking into account the previous discussion, the motion of an isolated ad- 
sorbate under the action of a bath consisting of the temperature-dependent 
surface vibrations can be characterized by the standard Langevin equation 

mR(t) = -m7R(t) + F(R{t)) + G{t). (6) 

Here, G(t) is the two-dimensional (Gaussian white noise) stochastic force, 
whose components satisfy the conditions given by equations ^ and ([5]) - 
moreover, {Gx{ti)Gy{t2)) = 0. Note that this equation of motion is based on 
the Markovian hypothesis: the correlations of the surface fluctuating force 
decay very rapidly [see equation ([5])] . On the other hand, also notice that 7 
gives rise to a characteristic timescale, the correlation time r = I/7, which 
is related to the mean free path that an adsorbate can travel without feeling 
much friction. 

In order to interpret the Langevin numerical simulations, we can express 
the intermediate scattering function, given by r.h.s. of the second equality 
of equation ([3]), as a second order cumulant expansion in A/C, 

I(AK, t) ^ e-^^' /,?(*-*')CAK(Odt' ^ (7) 

where C^i^it) = (wak(O) fAK(i)) is the autocorrelation function of the ve- 
locity projected onto the direction of the parallel momentum transfer. This 
is the so-called Gaussian approximation [12j, which is exact when the veloc- 
ity correlations at more than two different times are negligible. Despite its 
limitations, it provides much insight into the dynamical process by allowing 
an almost analytical treatment of the problem. Note that it allows to replace 
the average acting over the exponential function by an average acting over 
its argument, thus simplifying the analytical derivation in some simple but 
relevant cases. 

For an almost flat surface (T^ ~ 0), any direction is equivalent and there- 
fore the dimensionality of the numerical Langevin simulation reduces to one. 
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The corresponding numerical velocity autocorrelation function then follows 
the exponential behavior 

C{t) = (v^) e-^\ (8) 

where (uq) = ksTs/m (the square root of this magnitude gives the thermal 
velocity in one dimension). Introducing ([5]) into equation ([7]), we obtain [S] 

/(AK,t) =exp[-xMe-^* + 7t-l)] , (9) 
where x is the shape parameter, defined as 

X^^)^Kh = ~mK. (10) 



From this relation the mean free path results / = Ty^iv^, and the diffusion 
coefficient is D = t{v'^), which related to the friction (Einstein's relation). 
It can be easily shown [6J that the dynamic structure factor derived from 
the Fourier transform of equation ([9]) has a different shape depending on the 
value of X- This can be seen from the exact (analytical) Fourier transform 
of equation ([9]) , which renders the following functional form for the dynamic 
structure factor ruling the shape of the Q-peak 



5(AK,^) = x-'^Re|x-'*"^ r(x^ + i^r)-r(x' + i^r,x') 

(-l)V 2{x^ + n)/T 

n=0 



ex^^ (_l)n^2n 2{x^ + n)/T 

27r ^ n! uj^ + [{x'^ + n) / t]"^ ' ^ ' 



expressed in terms of the Gamma and incomplete Gamma functions (de- 
noted by r in the r.h.s. of the first line of equation (jlip to avoid any con- 
fusion with the FWHM, F), respectively. According to this expression, as x 
decreases (or, equivalently, as AK decreases, for Tg and 7 fixed) the shape 
of the dynamic structure factor goes from a purely Gaussian function to a 
Lorentzian one, and its width gets narrower and narrower (see appendix). 
This is the so-called motional narrowing effect, well-known in the theory of 
nuclear magnetic resonance lineshapes |6l[8l[9]. The shape parameter goes 
from zero to infinity. At high values of the Q-peak approaches a Gaussian 
shape, which will be well reproduced only by few terms (n not too small) in 
the sum given by equation (fTTI) . For values oix ^ 1, the Q-peak approaches 
a Lorentzian shape, with the term corresponding to n = usually being the 
dominant one in equation (jlll) . This variation between a Gaussian and a 
Lorentzian shape is in clear correspondence with having either a ballistic 
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Figure 1: Full width at half maximum, F, of the quasi-elastic peak as a 
function of the parallel momentum transfer, IS.K^ along the direction (100) 
of the Ft (111) surface for three different values of the friction coefficient: 
7 = (black/thick black line), 7 = 0.25 ps~^ (blue/thin black line), and 
7 = 2.0 ps^"*^ (red/thin grey line). 



(free particle) or a diffusive regime, respectively. Thus, a simple manner of 
expressing the FWHM of the Q-peak, F, in terms of x is 

F = 2^7X' + 2721^^2 (1-^)7X, (12) 

with ^ a free parameter. For /i = 0, one obtains the FWHM of the Gaus- 
sian lineshape, while the width corresponding to the Lorentzian lineshape is 
obtained for /x = 1 and when only the n = contribution of the infinite sum 
in equation (llip is considered (see appendix). That is, as the importance of 
the diffusive regime increases, one passes from a linear dependence on x to 
a quadratic one. The same behavior is observed in F when it is written as 
a function of IS.K because of the linear relationship between this magnitude 
and X (s6e equation p^ ). 

3 Results and discussion 

In figure [H F is plotted as a function of the parallel momentum transfer 
AK along the direction (100) for the Xe-Pt(lll) system at Ts = 105 K 
and three different values of the friction coefficient (0, 0.25, and 2 ps~^), 
assuming the corrugation of the Pt(lll) surface negligible. These results 
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Figure 2: Full width at half maximum, F, of the quasi-elastic peak as a 
function of the shape parameter for two different values of the friction co- 
efficient: 7 = 0.25 ps~^ (blue/thin black line) and 7 = 2.0 ps~^ (red/thin 
grey line). Vertical lines give the borders where a Lorentzian (L), a mixed 
Lorentzian-Gaussian (L-G), or a Gaussian (G) shape are observed. 



are in excellent agreement with those obtained by Ellis et al. [4| by means 
of Langevin molecular dynamics simulations, i.e. considering interacting ad- 
sorbates, though the interaction among them is relatively weak because the 
coverage used in the experiment was relatively low {6 = 0.017). The ex- 
perimental observation fitted perfectly the black solid line (see figure [l| , 
this confirming what they claimed as a two-dimensional free gas (zero fric- 
tion). In figure [1] it is also apparent a clear smooth transition in F, from a 
quadratic to a linear dependence on AK, in accordance with equation (fT2]) : 
the motional narrowing effect is clearly observed as the friction coefficient 
is increased for a given value of AX. The same behavior is also seen in 
figure O where F has been plotted as a function of the shape parameter. 
This figure also provides an important additional information: the gradual 
change of the shape of the Q-peak as a function of x- Very often, in order 
to extract information about the diffusion mechanism, a Lorentzian shape is 
assumed to deconvolute the experimental results. As clearly seen in figure[2l 
the pure Lorentzian shape is supposed to be good in a very narrow range of 
X (or IS.K) values. The three regions dividing this figure have been chosen 
according to the fittings of the numerical results obtained from the Langevin 
simulation to the analytical ones, given by equation (fTTIl . Though the two 
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perpendicular border lines denoting the different shapes are somewhat ar- 
bitrary, X values greater than one will give rise to a very strong Gaussian 
behavior. Therefore, this kind of representation has the advantage that ex- 
perimentalists can use it as a guide in the deconvolution procedure usually 
carried out within this context. 

The smallest experimental accessed value of the parallel momentum 
transfer, AET = 0.21 A~^, corresponds to a distance of about 30 A, which 
should be regarded as an upper limit to the mean free path of Xe atoms 
(or free gas regime) after reference [1]. Moreover, the experimental data for 
this system indicated that the Xe atoms run freely along the surface, with 
upper limits of 0.25 ps~^ for the friction and 9 meV for the diffusion barrier 
(at least less than the thermal energy). According to our theoretical study, 
the estimation of the mean free path should be carried out in a different 
way. The shape parameter and the friction can be easily obtained from a 
proper fit of the Q-peak to equation (jlip . and then the mean free path is 
finally computed using equation (jlOp . For example, for Ai^ = 0.21 
and 7 = 0.25 ps^^, one obtains that I = 3.2 A and that the correspond- 
ing lineshape of the Q-peak displays a mixture of Lorentzian and Gaussian 
shapes. 

The motional narrowing effect is also clearly seen in figures [3] and HI 
where numerical quasi-elastic lineshapes corresponding to different values of 
X (or Ai^) are shown. Dashed lines are fittings to the numerical lineshapes 
according to equation (jlip . In particular, in figure O the thick black curve 
corresponds to a zero friction case, for which the numerical (and analytical) 
lineshape is a pure Gaussian function. The blue curve is obtained from a 
Langevin calculation with a friction 7 = 0.25 ps^^, and the best fitting to 
equation pT]) is reached with rimax = 58, approaching a Gaussian shape. 
Conversely, the red curve displays a perfect Lorentzian behavior {n = 0, 
with a very slight correction of the lowest part of the Lorentzian that comes 
from n = 1). Analogously, in figure HI the numerical blue and red curves fit 
to a perfect Lorentzian (n = 0, with an also slight correction arising from 
71 = 1) and a mixed Lorentian-Gaussian function ijirnax 

= 10), respectively. 

4 Conclusions 

As seen in this work, the results obtained from the numerical Langevin sim- 
ulations carried out by Ellis et al. [4J fit perfectly to the analytical formalism 
based on the Gaussian approximation of the intermediate scattering func- 
tion. This allows to study the lineshape of the Q-peak in a low coverage 
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Figure 3: Dynamic structure factor, 5(AK,u;), as a function of uj for 
lS.K = 1.25 A^^ and three different values of the friction coefficient: 7 = 
(black/thick black line), 7 = 0.25 ps~^ (x = 4.1) (blue/thin black line), and 
7 = 2.0 ps~"^ (x = 0.5) (red/thin grey line). Dashed lines are fittings to the 
numerical lineshapes according to equation ([TT]) . 




CO (meV) 

Figure 4: Dynamic structure factor, 5(AK,u;), as a function of oj for 
7 = 0.25 ps"^ and two different values of AA': AK = 0.15 A"^ (x = 0.5) 
(blue/thin black line) and Ai^ = 0.5 A^^ (x = 1-6) (red/thin grey line). 
Dashed lines are fittings to the numerical lineshapes according to equa- 
tion (HI]). 
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regime and for very low corrugated surfaces, as it is the case for the Pt(lll) 
surface. Accordingly, we have shown that a gradual change of the shape 
of the Q-peak is expected when increasing the friction and AX, which can 
be understood in the light of the so-called motional narrowing effect. This 
type of analysis of the Q-peak should be therefore prior to any deconvolu- 
tion procedure carried out in order to better extract information from the 
process (e.g. diffusion coefficients and jump distributions). Moreover, taking 
into account these numerical simulations, conditions for the corresponding 
free gas regime quite different to those reported earlier have been found. 
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Appendix 

In this appendix we analyze the two extreme cases of the dynamic struc- 
ture factor lineshape (Lorentzian and Gaussian) as a function of the shape 
parameter, starting from equation (jllh . 

For X 1) the dominant term in the second equality of equation (jlip is 
that corresponding to n = 0. Therefore, one obtains straightforwardly that 

SiAK,u)^-^ \ (13) 

which is a Lorentzian function with FWHM F = 27%^ = 2DAK^ (i.e. fJ- = 1 
in equation p^ ). 

On the other hand, for large Xi it is convenient to start with the first 
equality of equation (|11|) , which is written in terms of the incomplete Gamma 
and Gamma functions [13j, whose asymptotic behaviors are 

f(a,/3) ^ P'^-'e-^, (14) 
f(Q) « V2^a"-i/2g-a^ (15) 
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respectively, with a = + ^"^/t and /? = • For the sake of simphcity, it 
is better to express a is its polar form, a = pe^^, with 



ijj 
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b = (tan)' 



1 ( 



(16) 



(the approximation in b arises after assuming large x). In doing so, equa- 
tions ([HI) and (fT5l) become 



f(x2 + i../7,x')-X^(^^+^'^/^)e-^\-2 



(17) 



and 



V{x + i^h) = V27r 



x'.(^)- 



1/2- 



+ 1 - 

7 



-1/4 



(18) 



respectively. After some straightforward algebraic manipulations, equa- 
tion (I18D can be rewritten as 



(19) 



Substituting (fT7|) and (fT9|) in equation pT]) we finally obtain 

-1 -a;V272x2 



5(AK,( 



1 r 



7 



X 



which in the limit of large x reads as 



5(AK,cj) 



27r 



7X 



-a;2/27'x' 



(20) 



(21) 



This is a Gaussian function with F = 2^/2 In 2 7x esc Ai^' (i.e. n — > oo in 
equation pT|) and /i = in equation (fT2|) ). 
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